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We propose a new platform for interacting topological phases of fermions with time reversal 
symmetry © (such that ©^ = 1) that can be realized in vortex lattices in the surface state of 
a topological insulator. The constituent particles are Majorana fermions bound to vortices and 
antivortices of such a lattice. We explain how the © symmetry arises and discuss ways in which 
interactions can be experimentally tuned and detected. We show how these features can be exploited 
to realize a class of interaction-enabled crystalline topological phases that have no analog in weakly 
interacting systems. 


Theoretical understanding of topological phases of 
non-interacting fermions, now thought to be complete, 
gave us a treasure trove of new materials over the past 
decade including the topological insulators (both 2D and 
3D) [IHl] and engineered topological superconductors 
with Majorana zero modes [SHE]. The same theoreti¬ 
cal framework also opened new areas of inquiry, leading 
to concepts such as the Weyl semimetals, axion insula¬ 
tors and topological exciton condensates that are actively 
pursued in experiments. When strong interactions are in¬ 
cluded the situation becomes more complicated and the 
classification of topological phases less well understood, 
especially for fermionic systems [9]. It is expected that 
interacting systems could produce novel phases with un¬ 
usual properties, such as excitations with fractionalized 
quantum numbers and unusual exchange statistics. Al¬ 
though a number of specific interaction-enabled topolog¬ 
ical phases with fermions have been theoretically pro¬ 
posed [IDHIH], the only such phases that are known to 
exist in the real world are the fractional quantum Hall 
states [191 [20]. 

With the goal of enlarging the space of experimentally 
accessible topological phases that depend for their exis¬ 
tence on strong interactions we propose here a physical 
platform that realizes the paradigm of fermions with time 
reversal (TR) symmetry 0 such that 0^ = 1. Interact¬ 
ing fermions with this property have been employed in 
seminal works by Fidkowski and Kitaev (FK) [21] who 
showed how the integer classification of their ID topo¬ 
logical phases in the non-interacting limit changes to 
Zg classification when interactions are included. More 
recently, Lapa, Teo and Hughes (LTH) [22] introduced 
a model with 0 and an additional inversion symmetry 
V. In this case, it turns out that there are no topologi¬ 
cally non-trivial ID phases in the absence of interactions. 
Remarkably, when interactions are included, a topologi¬ 
cally non-trivial phase becomes possible (the interacting 
classification here is Z 2 ). The latter, then, is a genuine 
interaction-enabled topological phase that fundamentally 
cannot exist in the non-interacting limit. 

The TR symmetry that we have in mind acts on the 


spinless fermion annihilation operator Cj through an op¬ 
erator 0 so that 0Cj0“^ = Cj and 0i0“^ = —i. If we 
decompose our Dirac fermion into a pair of Majorana 
fermions Cj = \{oLj + i^j) it follows that 

QajQ~^ = aj, = —^j. ( 1 ) 

In a ©-invariant Hamiltonian expressed in terms of aj , pj 
operators only certain terms are allowed. For instance bi- 
linears iajPk ^-re allowed while ia^a^ and iPjPk ^-re pro¬ 
hibited. Similarly, four-fermion interaction terms with 
even numbers of ajA are allowed, such as ajak^iam or 
PjPkC^iam^ but those with an odd number are not. Find¬ 
ing a physical system whose Hamiltonian implements the 
above symmetry constraints presents a challenge because 
electrons, the only relevant fermions in solids, have spin ^ 
and the natural TR operation 0 for such spinfull fermions 
satisfies 0^ = —1. The important ideas [2TH23] that in¬ 
volve fermions with 0 such that 0^ = 1 have therefore 
remained largely untested (see however proposal in Ref. 
[24]). In what follows we show how the ingredients nec¬ 
essary to realize the phases envisioned by FK and LTH 
can be realized in a concrete physical system accessible 
with the available experimental methods. 

A system we explore in this work is similar to that 
studied in Ref. [25] - a superconducting (SC) surface of 
a strong topological insulator (STI) - with one extra in¬ 
gredient. In addition to vortices, which are known to 
harbor unpaired Majorana zero modes (MZMs) [26], we 
include in our considerations antivortices which also con¬ 
tain MZMs but, as we show, of a different type. Specif¬ 
ically, we demonstrate that it is consistent to assign the 
two types of Majorana fermions aj, pj that obey Eq. 
0 to vortices and antivortices, respectively. Structures 
composed of vortices and antivortices, arranged such 
that their corresponding Majorana wavefunctions have 
non-zero overlaps, then implement ©-invariant fermionic 
Hamiltonians with 0^ = 1. We discuss various vor- 
tex/antivortex geometries that realize interacting lattice 
models, including the LTH model [22] mentioned above. 

We now proceed to substantiate these ideas and claims. 
The physical system, an STI with a superconducting sur- 
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face, is described by the Fu-Kane Hamiltonian [26] 1-L = 
/ (fr^lHFK{r)^r, where 4 '^ = (c-|-^, 04,^, c|^,is 
the Nambu spinor and 

^FK = r^{p • tr - /i) + r^Ai + r^A2. (2) 

Here tr, r are Pauli matrices in spin and Nambu spaces, 
respectively, and A = Ai + i/S .2 represents the SC or¬ 
der parameter. A single isolated vortex, expressed as 
A(r) = Ao(r)e“*’^^ with (p the polar angle and n = ±1 
corresponding to vortex or antivortex, respectively, is 
known to bind a MZM [26] . The Hamiltonian Q respects 
the particle-hole symmetry C generated by S = a^r^K 
= -1-1, K denotes complex conjugation) and, for a 
purely real gap function A, also the physical TR sym¬ 
metry 0 generated by 0 = ia^K (0^ = — 1). In the 
presence of vortices 0 is broken but in the special case 
when /i = 0, the Hamiltonian respects a fictitious TR 
symmetry 0 with 0 = a^r^K (0^ = +1), even in the 
presence of vortices m- From now on we focus on this 
/i = 0 “neutrality point” with an extra symmetry. To¬ 
gether, the two symmetries E and 0 define a BDI class 
with chiral symmetry H = S0 = 

Eigenstates of the Hamiltonian © are four-component 
Nambu spinors ^(r) = Because the 

particle hole symmetry C maps positive energy eigen¬ 
states to their negative energy partners, a non-degenerate 
zero mode H^k^o = 0 is self-conjugate under C, that 
is, it obeys = ^o- This constrains its components 
such that and = —u^. In addition, because 

{i^FK^n} = 0, the zero mode must be an eigenstate 
of H. There are two choices that satisfy these conditions, 

= = (3) 

corresponding to the two eigenvalues z/ = ±1 of H. It is 
easy to check that 


(4) 

i.e. the two MZMs transform as even and odd under 0. 
To complete the argument we construct the correspond¬ 
ing zero mode operators j± = f From 

Eqs. (© and @ it follows that 

7± =7±, 07±0“^ = ±7±. (5) 

The zero mode operators are Majorana and the two types 
transform in the opposite way under 0. Now suppose 
that 7 + with a wavefunction resides in the core of 
a positive n = 1 vortex (as can be verified by an ex¬ 
plicit calculation [25]). The physical TR symmetry 0 
maps a vortex onto an antivortex because the direction 
of superflow is reversed under 0. At the same time the 
corresponding operator 0 maps to \ Antivortex 
thus necessarily carries the other type of Majorana rep¬ 
resented by 7 _. Comparing Eqs. (§ and 0 we conclude 
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FIG. 1: Vortex lattice geometries for interacting Majorana 
models. Solid (open) circles represent vortices (antivortices), 
while solid and dashed lines indicate hopping and interaction 
terms, respectively, consistent with ©. a) Cluster of 8 vortices, 
b) The 4-leg LTH ladder invariant under both © and V. Pos¬ 
sible experimental realizations are sketched in panels c) and 
d). In panel (c) the STI surface is assumed to be supercon¬ 
ducting and the thickness of the STI sufficiently large so that 
the Majorana wavefunction overlaps occur predominantly in 
the surface. Under these conditions the setup will realize the 
LTH ladder depicted in Fig.[^. In panel (d) the interactions 
take place in the surfaces but tunneling is assumed to occur 
through the bulk of the flake. 


that vortices and anti vortices in the Fu-Kane model at 
neutrality carry MZMs that transform as even and odd, 
respectively, under 0 and can thus be assigned as aj , Pj 
Majorana operators in models with 0^ = I. We note 
that these results can be explicitly checked by a direct 
calculation of the zero-mode wavefunctions [8] in Hamil¬ 
tonian 0. 

Systems in this symmetry class composed of only vor¬ 
tices and no antivortices admit interaction terms of the 
form gajakO^iam but no bilinear terms and are thus in¬ 
herently strongly interacting. Some models of this type 
have been explored in Ref. [25]. When the system is 
slightly detuned from neutrality then 0 is broken and 
bilinears it'aja^ become allowed with t' ^ ji. In the 
following we shall adopt an assumption that p has been 
tuned sufficiently close to zero so that t' g and we may 
thus neglect all bilinears prohibited by 0. 

As our first example of interesting structures that can 
be constructed with these ingredients we consider a clus¬ 
ter made with a small number n of vortices. As we shall 
argue, scanning tunneling microscopy (STM) can be used 
to probe the effects of interactions in such a small cluster 
by examining the vortex core spectra for the presence or 
absence of MZMs. We envision changing the total vortic- 
ity n of our cluster from 0 to 8 by adding vortices one by 
one to observe the theoretically predicted Zg periodicity 
[2T1 [23] generated in the presence of interactions. In the 
absence of interactions the ground state has degeneracy 
2^^/^ and STM will observe a single MZM in each vortex. 
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FIG. 2: Ground state degeneracy in a cluster with n vor¬ 
tices and generic interaction terms. Stars denote presence or 
absence of zero bias peaks in vortices observable by STM. 


When interactions are present, STM should still see zero 
modes for n = 1, 2,3 and for all odd values of n, but the 
zero modes could generically split for even n > 4, be¬ 
cause an interaction term can be first constructed with 4 
Majorana operators. The pattern of ground state degen¬ 
eracies we obtain for a system of n vortices with generic 
interactions allowed by 0 is displayed in Fig. and is 
consistent with results of Ref. We find that MZMs, 
detectable by STM, will be present for all n except when 
n = Ak with k integer. This is to be contrasted with a 
non-interacting case with generic hopping terms allowed 
(e.g. when /i ^ 0); here the ground state degeneracy is 
^/ 2 {n mod 2) and zero bias peaks will be seen for all odd 
n. 


As an example of considerations that lead to Fig. 
we now discuss the special cases of an interacting system 
with n = Ak. We start with n = 4, denote the Majorana 
operators as o^i, 0 ^ 2 , 0^2 5 define complex fermions 
di = ^{ai A- ia'i). Noting that iaia[ = 2ni — 1 where 
ni = djdi is the number operator we may label the quan¬ 
tum states of the 4 Majorana fermions in the cluster by 
the eigenvalues ni =0,1 as |nin 2 ). The most general 
interaction term = gaia 2 a'iCX 2 splits the 4-fold de¬ 
generacy into an even parity ground state doublet |00), 

111) (for ^ > 0) and odd parity excited states |01), |10). 
Because both ground states have the same parity single¬ 
electron tunneling will necessarily cause transitions to the 
excited states and the STM peaks will appear at energies 
Az 2 g^ not zero. 

For n = 8 we consider a specific pattern displayed in 
Fig. This pattern is easy to analyze and also forms 
the basic building block for the interaction-enabled LTH 
topological crystalline phase that we shall discuss below. 
As argued previously [25] strongest interactions occur for 
those groups of 4 Majorana fermions that are spaced most 
closely together. This is because the corresponding cou¬ 
pling constants depend on the overlap of the exponen¬ 
tially decaying Majorana wavefunctions. With this in 
mind we identify the dominant interaction terms associ¬ 


ated with the square plaquettes in Fig. described by 

h\j = gi{aiaia2Q^2 + a3a'^a4a'^) + ^2 <^2 <^2 ^3 <^ 3 - ( 6 ) 

For ^ 1,^2 > 0 the ground state of h\j is doubly degen¬ 
erate, spanned by eigenvectors | 0000 ) and in the 

same notation as above. The energy is Eg = —2gi — g 2 . 
The sub dominant interaction term is associated with two 
straight legs. 


h\ = g'{aia 2 a 3 a 4 + (7) 

It is easy to see that the inclusion of h\ splits the ground 
state doublet of h\j into a bonding/antibonding pair with 
= (|0000) ± \llll))/^/2. For g' > 0 the unique 
ground state of 


^8 = + h\ ( 8 ) 

is IV^-), while |' 0 +) is the first excited state with energy 

V. 

Now consider the 4-leg chain depicted in Fig.[^. This 
is a version of the LTH model [22] that can be constructed 
using our platform. It consists of alternating clusters of 
8 vortices and antivortices, each described by interact¬ 
ing Hamiltonian hs{cx 2 j) and hs{l 32 j-\-i), connected by 
nearest-neighbor hopping terms with amplitude t. Here 
OL 2 j denotes the octet of {ai^a[} operators in cluster 2 j 
and the same for I32j+i- The Hamiltonian, then, reads 

N N 

^LTH = ^ ^ ^8(/32i+l) + 

j=-N j=-N-l 

AT+l 

- it ^ { 0 i'l^ 2 jPl, 2 j+l + Pl 2 j-l<^l, 2 j), (9) 

j=-N-l 

where has a form indicated in Eq. (j^ and describes 
the four dangling Majoranas at the two ends of the chain. 
The Hamiltonian ^ respects 0 as well as the inversion 
symmetry V. The latter is generated by {aij^a'ij} 

{a[ _j,ai_j} and A'j} ^ -{/3; 

Without interactions, there exist no topologically non¬ 
trivial phases in a system with these symmetries [ 22 ]. 
Indeed we see that when gi = g 2 = 9 ' = ^ and t 7 ^ 0 the 
system breaks up into a set of local dimers formed by 
products on the dashed bonds in Fig. [^. The ground 
state is unique with a gap 2 t to the lowest excitation and 
clearly topologically trivial. 

In the opposite limit, t = 0 and gi g 2 > g' > 0, the 
ground state is a direct product of \ip-)j states on each 
cluster j. If we impose periodic boundary conditions the 
ground state is unique with a gap Ag' to the lowest ex¬ 
cited state. However, for an open chain that preserves 
the symmetries, as the one indicated in Fig. [^, there is 
a 4-fold degeneracy associated with The four in¬ 

dependent quantum states associated with the quartet of 
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FIG. 3: Excitation energies of the LTH Hamiltonian § ob¬ 
tained by exact numerical diagonalization. a) The lowest 
excitation energy in a system of length Lx = 4,8,12 with 
periodic boundary conditions, b) Several lowest excitation 
energies for = 8 and open boundary conditions. We take 
9i — 92 — g' = g and the energies are in the units of t. 

ai Majoranas at each end are split by into a dou¬ 

bly degenerate ground state (with even local parity) and 
a doubly degenerate excited state (with odd parity). As 
noted in Ref. [22] , 0 connects the two degenerate ground 
states, but its action in this subspace is anomalous with 
0^ = —1. The edge modes therefore comprise two effec¬ 
tive spin-^ degrees of freedom and constitute fractional¬ 
ized excitations analogous to those appearing in spin-1 
Haldane chains [28l [29] . In a long chain the edge degen¬ 
eracy is therefore protected by the Kramers theorem and 
cannot be removed by any local perturbation preserving 
0. It signals a topological phase, one that fundamentally 
cannot exist in a non-interacting system. If we turn on 
a small hopping t the topological phase should persist 
up to a critical strength tc of order g' ^ at which point a 
phase transition occurs to the trivial phase. 

We have confirmed the above picture by performing 
exact numerical diagonalizations of F^lth; the results 
are displayed in Fig. Our simulations with periodic 
boundary conditions (panel a) indicate a phase transition 
marked by the excitation gap closing at Qc — 0.81t, ob¬ 
tained by extrapolating the energy minimum to ^ oo. 
Panel b confirms that with open boundary conditions 
the ground state is unique for g < gc but becomes 4- 
fold degenerate for g > gc with all ground states in the 
same parity sector, in accord with our expectation for 
the interaction-enabled topological phase. 

As explained in the Supplementary Material m a sys¬ 
tem described by the LTH Hamiltonian ^ can be re¬ 
garded as a “4e superconductor”. (This is because it 
breaks the fermion number conservation symmetry while 
all the anomalous expectation values that are bilinear, 
such as (di(i 2 ), vanish.) The authors of Ref. [22] proposed 
to look for physical realizations of such a 4e superconduc¬ 
tor in certain pair density wave systems, discussed the¬ 
oretically in the context of underdoped high-Tc cuprate 
superconductors [31]. Whether or not such a phase ex¬ 
ists in real materials remains to be seen. Our proposed 
realization of the LTH model, on the other hand, relies 


on ingredients that are known to exist. Specifically, SC 
order in the surface state of an STI has been experimen¬ 
tally observed by a number of groups [32^4^ . In some 
cases the chemical potential has been tuned to the close 
vicinity of the Dirac point [39] [40] as required for models 
with symmetry 0. Vortex cores m, and the expected 
MZMs have also been imaged [43]. Assembling vortices 
into regular structures, controlling their geometry and 
reliably probing the zero modes remains a challenge but 
it is one that does not seem insurmountable given the 
recent progress. 

In the near term it should be possible to probe the 
effects of interactions in small clusters of vortices by 
STM, as we discussed. The relevant energy scale can 
be quite large, ^ lOmeV, under favorable conditions [25] 
and should permit observation of the interaction induced 
Zg periodicity. To engineer the LTH model the key chal¬ 
lenge will be to assemble stable arrays containing both 
vortices and antivortices in close proximity to one an¬ 
other. We note that such structures have been observed 
to occur spontaneously in mesoscopic SC samples with 
certain geometries [44] . Alternately, one can leverage the 
fact that in a thin STI/SC film or flake a uniform per¬ 
pendicular field B will produce vortices on the top sur¬ 
face and antivortices on the bottom surface. Panels (c) 
and (d) in Fig. [^outline two possible experimental real¬ 
izations of the LTH vortex/antivortex lattice exploiting 
this principle [45] . 

The interaction-enabled topological phases discussed 
in this work are gapped and therefore robust with respect 
to moderate amounts of symmetry-preserving disorder. 
Specifically, disorder in vortex positions does not break 
0 and is therefore innocuous. It also follows that fine de¬ 
tails of the geometry sketched in Fig.[^a,b) are unimpor¬ 
tant as long as 0 and V are preserved to a good approx¬ 
imation. Any experimental realization that reasonably 
approximates the proposed model geometry should show 
the topological phase. Local fluctuations in the chemi¬ 
cal potential break 0 but we expect the gapped phases 
to remain robust as long as the symmetry is preserved 
on average and the fluctuations do not exceed the gap 
amplitude HIIB]. 
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SUPPLEMETARY MATERIAL 


Disorder effects 


As mentioned in the main text we expect symmetry¬ 
preserving disorder to have only a mild effect on the 
gapped topological phase of the LTH model. Here we 
present numerical evidence that supports this expecta¬ 
tion. We note that the roles played by the two sym¬ 
metries 0 and V are somewhat different. The former 
underlies the existence of the topological phase and pro¬ 
tects the exact degeneracy of its edge modes. The latter 
plays a subsidiary role by prohibiting the existence of the 
topological phase in the absence of interactions. Con¬ 
sequently, we expect ©-breaking disorder to have more 
pronounced effects on the system than the 7^-breaking 
disorder. This is indeed borne out by numerical simula¬ 
tions presented below. 

To simulate the 7^-breaking disorder we replace the 
hopping t ^ and the interaction g ^ g ^ 5g in the 

LTH Hamiltonian Eq. (9). Here 5t is a random variable 
drawn independently on each bond from a uniform distri¬ 
bution 5t G {—Wt^Wt) and similarly ioi 5g ^ {—Wg^Wg). 
Fig. 4 shows the lowest excitation energies of the disor¬ 
dered system with open boundary conditions. We ob¬ 
serve that the gap and the expected 4-fold ground state 
degeneracy in the topological phase remain robust for 
moderate levels of 7^-breaking disorder (top panels in Fig. 
4). To simulate disorder that breaks 0 we add 6 H to the 
LTH Hamiltonian Eq. (9) that contains terms i5t'aija[ j 
and i6t'j describing tunneling events between Ma- 
jorana sites within each cluster. Such terms are prohib¬ 
ited by 0. The amplitude is again drawn from a random 
distribution 6 t' G The results, displayed in 

bottom panels of Fig. 4, show that while the gap remains 
robust with respect to moderate levels of ©-breaking dis¬ 
order the exact degeneracy of the edge states is now lifted 
in proportion to the disorder strength. This is akin to the 
gap formation observed in the surface state of a topologi¬ 
cal insulator induced by doping with magnetic impurities. 
Importantly, the edge modes remain nearly degenerate 
and clearly distinguishable from the bulk modes even for 
a fairly large disorder strength w[ = O.lt. 

We conclude that for all practical purposes the 
interaction-enabled topological phase survives the inclu¬ 
sion of moderate levels of both 0- and 7^-breaking disor¬ 
der as expected on the basis of general arguments. The 
bulk gap and the gapless edge modes remain experimen¬ 
tally observable in the presence of disorder. In fact these 
signatures of the phase should survive until the disorder 
strength becomes comparable to the bulk gap. 



FIG. 4: Excitation energies of the LTH Hamiltonian with 
disorder obtained by exact numerical diagonalization in a sys¬ 
tem of length Lee = 8 with open boundary conditions. The 
same model parameters as in Fig. (3) are taken. The top 
row shows results for 7^-breaking disorder with strength wt = 
Wg = 0.05t and Wt — Wg = O.lOt. The bottom row shows 
results for ©-breaking disorder with strength w't = 0.05t and 
w't = O.lOt. Because these simulations are numerically costly 
we do not average over disorder realizations. However, each 
point in these graphs corresponds to an independent disor¬ 
der realization which provides a clear sense for the statistical 
properties of the system. 


Comments on the LTH model realization 

Fig. 1(b) represents an ideal configuration for the real¬ 
ization of the physics embodied by the LTH model. Some 
experimental geometries that can be used to approximate 
this configuration are sketched in panels (c,d) of that Fig¬ 
ure. In these, red (blue) dots represent vortices (antivor¬ 
tices) and it is assumed that the lateral extent of both 
structures is such that it accomodates 4 vortices and an¬ 
tivortices in the direction perpendicular to the plane of 
the drawing. 

Aside from the dominant couplings {gi^g 2 ^g' ^ t) consid¬ 
ered in the model Hamiltonian Eq. (9) in the real physi¬ 
cal system all other couplings allowed by symmetry will 
also be present. For instance there will be additional 
hoppings across the plaquette diagonal as well as addi¬ 
tional four-fermion terms within each cluster and those 
involving two a and two p operators between the clusters. 
Although not explicitly included in our numerical calcu¬ 
lations, such symmetry-preserving terms will not alter 
the conclusions regarding the interaction-enabled topo¬ 
logical phase in this model, as long as they remain small 
and short ranged. The latter property is guaranteed by 
the exponential decay of MZM wavefunctions while the 
former must be built into the geometry of the structure. 
The key point here is that a generic set of four-fermion 
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FIG. 5: Excitation energies of the LTH Hamiltonian as a 
function oi g — gi — g 2 — 3g'. 

interaction terms acting within a cluster of 8 MZMs pro¬ 
duces a non-degenerate ground state, except for acciden¬ 
tal degeneracies that occur on a set of measure zero in 
its parameter space. This is, ultimately, all we need to 
guarantee the existence of the gapped topological phase 
in the limit of strong interactions. Inclusion of the sym¬ 
metry preserving hopping terms (again, of an arbitrary 
form) will then sustain the gapped phase for some range 
of hopping strength as illustrated in Figs. 3 and 4. To 
further document this robustness we have performed nu¬ 
merical computations using the LTH Hamiltonian (9) for 
different sets of model parameters. As an example we 
show in Fig. 5 the energy spectrum of the system for the 
case when g' is smaller than both gi and ^ 2 - As expected 
on general grounds this leads to the same qualitative re¬ 
sults for the phase diagram of the system. We also note 
that, for the model to describe spinless fermions with 
time reversal symmetry, we require equal numbers of a- 
and /3-type MZMs as well as certain notion of locality, i.e. 
that MZM components of the complex fermion Cj defined 
above Eq. (1) be not too far removed from one another. 

These observations, combined with the robustness 
against disorder demonstrated above, indicate consider¬ 
able leeway in the actual experimetal implementation. 
The experimental system need not have the precise geom¬ 
etry indicated in Fig. 1(b); so long as it consists of alter¬ 
nating clusters of 8 vortices and antivortices, connected 
with one another predominantly by tunneling terms, and 
so long as it preserves the V and 0 symmetries to a 
reasonable approximation, the interaction-enabled topo¬ 
logical phase should naturally obtain in the regime when 
interactions dominate over the tunneling terms. 

While the geometry is rather flexible there are some ba¬ 
sic requirements on the parameters of the physical system 
that must be met for it to realize the interacting phase. 
Most importantly, we require that |/i| is sufficiently close 
to zero so that the direct tunneling t' ^ ji between MZMs 
of the same type is negligible compared to the interac¬ 
tion scale g. We have previously estimated g lOmeV 
under favorable conditions which means that /r must be 
tuned to zero to within few meV. Second, we require that 
MZMs are well separated in energy from all other vortex 
core states. The lowest lying such state has the energy 


El ^ /S? / -f- jE. We see that when /i C A this en¬ 

ergy approaches the full SC gap so the second condition 
ultimately also boils down to the chemical potential be¬ 
ing as small as possible while simultaneously maintaining 
significant gap amplitude A. 

LTH model as an emergent 4e superconductor 

As mentioned in the main text the gapped topological 
phase of the LTH model can be regarded as a 4e super¬ 
conductor (SC). Here we explain this notion in greater 
detail. First, the system should be considered a super¬ 
conductor because its Hamiltonian (9) breaks the parti¬ 
cle number conservation symmetry (although it conserves 
the fermion parity). In an ordinary SC such number non¬ 
conservation is associated with a non-vanishing anoma¬ 
lous expectation value = {d\d\,) often identified as 
the SC order parameter. Here, di = Eici'i) is a com¬ 
plex fermion operator defined in the main text. However, 
deep in the topologial phase of the LTH model all such 
anomalous expectation values vanish. This can be seen 
most easily by considering the limit gi = g 2 ^ g' and 
t = 0 of the model. As discussed below Eq. (7) the 
ground state in this limit can be represented as a direct 
product 

l^> = n,-|V’-)i, (10) 

where = (|0000)j ± |llll)j)/v^ are the two low 

lying states of the 8-site cluster j. The low-lying excited 
states are obtained by replacing some in Eq. ( pT] ) 

by |'0+)j and have energy Ag' for each excited cluster. 

Because the two basis states |0000) and differ 

in occupancy by 4 fermions it is clear that = 0 

for all /,/' when evaluated in the ground state |T). It 
is also clear that a 4-fermion expectation value A^^^ = 
{d\d\d\d}^ will be non-zero (of order one) in this ground 
state. Away from the t = 0 limit the ground state does 
not have a simple form but, at least for small t, can be 
written approximately as a linear combination of prod¬ 
uct states constructed from As long as this is the 

case the property = 0 and A*^^^ = C)(l) continues to 
hold and the system must be regarded as a 4e supercon¬ 
ductor. Only once t grows comparable to g' the ground 
state begins to contain a small admixture of higher en¬ 
ergy cluster states, such as |0011), possibly causing small 
non-zero values of . Nevertheless, the essential char¬ 
acter of the topological phase remains encoded in A*^^^ 
all the way to the phase boundary with the trivial phase. 








